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Abstract 

In this article, we study the self-similar solutions of the 2-component Degasperis-Procesi 
water system: 

pt + k-zupx + (fci + ki)pu x = 
lit — Uxxt + 4uu x — 3u x u xx — uu xxx + kzpp x = 0. 
By the separation method, we can obtain a class of self-similar solutions, 



(1) 



p(t,x) = ma X ( a(4t)( { ( ^ fc2)/4 , 0), u{t,x) = 
a ( s ) ~ 3^ivF = °> a (°) =00^0, a(0) = ai 



o(4t), 



4a(s) K 



(2) 



where r\ = ^ 1/4 with s = 4i; k = + &2 — 1, a > 0, £ < 0, ao and ai are constants, 
which the local or global behavior can be determined by the corresponding Emden equation. 
The results are very similar to the one obtained for the 2-component Camassa-Holm equa- 
tions. Our analytical solutions could provide concrete examples for testing the validation and 
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stabilities of numerical methods for the systems. With the characteristic line method, blowup 
phenomenon for £3 > is also studied. 
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1 Introduction 

The 2-component Degasperis-Procesi shallow water system [20] and [10], can be expressed in the 
following form 

p t + k 2 up x + (ki + k 2 )pu x = 0, x e R, 

(3) 

ut - u xxt + iuu x - 3u x u xx - uu xxx + k 3 pp x = 0. 
Here k±, k 2l k^ are constants, u = u{x,t) € R is the velocity of fluid. And p — p(t,x) > is the 

density of fluid. For p — 0, the system returns to the Degasperis-Procesi equation [7], [17] and [26] . 

In this article, we adopt an alternative approach (method of separation) to study some self- 
similar solutions of 2-component Camassa-Holm equations ([3]). Indeed, we observe that the isen- 
tropic Euler, Euler-Poisson, Navier-Stokes and Navier-Stokes-Poisson systems are written by: 

Pt+V ■ (pu) =0 

( p[u t + (w • V) u)] + VP = — pV$ + vis(p, u) ( 4 ) 
A$(t,x) = a(N)p 

where a(N) is a constant related to the unit ball in R N : a(l) — 2; a(2) — 2tt and for N > 3, 

a{N) = N(N - 2)Vol{N) = N(N - 2) + (5) 

where Vol(N) is the volume of the unit ball in R N and T is a Gamma function. And as usual, 
p = p(t, x) and u = u(t, x) £ are the density and the velocity respectively. P = P{p) = Kp r 
is the pressure, the constant K > and 7 > 1. And vis(p,u) is the viscosity function. 
We may seek the radial solutions 

p(t, x) = p(t, r) and u = -V(t, r) =: -V (6) 
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with r = (J2? =1 , 

By the standard computation, the Euler equations in radial symmetry can be written in the 
following form: 

f AT — 1 

I Pt + Vp r + pV r H pV = 

(7) 



p (V t + vv r ) + KJLp-y = 0. 

For the mass equation in radial symmetry, ([T]h . we well know the solutions' structure (Lemma 3, 

my- 

f(-^) nM 

(8) 



, . /vSTTf) / \ a(r) 

p^ r ) = -mr> < l > r )=-K r - 



As the 2-component Degasperis-Procesi equations ([3]), are very similar to the Euler system ([7]), 
we can apply the separation method ([13], [IB], [H], [21]) to the systems ([SJ- We note that 
for the 2-component Camassa-Holm equations [Tl], [I], [T2] and [5], 



pt + up x + pu x = 
m t + 2um x + urn^ + up/Or = 



(9) 



with 



m = ii — a Uj 



(10) 



By the separation method, we can obtain a class of blowup or global solutions for a = 1 or — 1 
[25] . In particular, for the integrable system with a = 1, we have the global solutions: 



KM) 



°(3t) . 
a(3t)" 



<3t)T73> for »T < 
0, for ?y 2 > £a 2 

a ( s ) ~ 3a( f)i/3 = °> a (°) = a o > 0, a(0) = a i 



(11) 



m = |\/-^ 2 +(ea) 2 

where r\ = a ^ij 3 with s = 3t; £ > and a > are arbitrary constants. 

As the two component of Camassa-Holm equations are very similar to the two component of 
Degasperis-Procesi equations, In this article, we continue to apply the separation method, to 
deduce the nonlinear partial differential equations into much simpler ordinary differential equations. 
Therefore, we can contribute a new class of self-similar solutions in the following corresponding 
result: 



4 



M.W. Yuen 



Theorem 1 We define the function a(s) is the solution of the Emden equation: 

a(s) - t4vt = 



4-(a)« (12) 



a(0) = a > 0, a(0) = a 



and 



«">-TV-5 , " + (f^' (13) 



where r\ — a ^ k2 /i with s = At; K = 4r + k 2 — 1, a > 0, £ ^ 0, ao and ai are constants. 

For the 2-component Degasperis-Procesi equations {5}), i/iere exists a family of solutions, those are: 

(l)for ks = 0, and £ = 0, or 

P ( tx ) = ^ uftx) = ^ a: (14) 



where po > 0, is an arbitrary C 1 function. 

(2) for fc 3 > and f > 0, or 

(3) for fc 3 < and £ < 0, 



P ft >g )=max( o(4 ff + ) fca)/4 , 0) , u(t,x) = fjgg*. (15) 



Remark 2 TTie structure of solutions U5\) of the 2-component Degasperis-Procesi equations are 
very similar the one ill}) . [25] of the Camassa-Holm equations. 

2 Separation Method 

The mass of the solution is not conserved except on k\ = and k 2 — 1 in equation ([3j) i . However, 
we can also design a nice functional structure for the mass equation: 

Lemma 3 For the l-dimensional equation of mass f2pi.' 

pt + k 2 up x + (ki + k 2 )pu x = (16) 

there exist solutions, 

i± \ ^aut) k 2/i) . a(At) 

p(t,x)= , .,, ' , w , , u(t,x) =^7ttx (17) 

with the form f(r\) > G C 1 wii/i n = a(4t ^ 2/4 , and a(4i) > G C 1 . 
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Proof. We just plug (JT7J) into dTSJ) to check: 



p t + k 2 up x + (ki + k 2 )pu x 
d I /( 0(4*^2/4) \ 



(18) 



, a(4t) d ( /( or^a/*) 
fc 2 , ' 



St I a(4i)( fe i+ fe 2)/ 4 / a(4i) & 1 o(4f)( fc i+^)/ 4 



(fci + fc 2 ) 



f( a (4t) k ''*) d ( a(4t) 



!(4t)(fci+fca)/4 ^ I a(4*)* 



(19) 



a ( 4i )(fel+fc 2 )/4- 



(fci + fc 2 ) . Z 1 • X _g X 

- 1 1 4 ' ' a ^ ' 4 ' /l a (4<)^/4 J + a ( 4 <)(fei+fe2)/4 / \(4t) k ^ 1 9r a (4t) fe 2/4 ' 



fc 2 a(4i):r /( a( 4t)*2/4 ) 9 ^ a: ^ ( fc i + M/( a(4t )fc 2 /4 ) q(4£) 

a(4t)( fe i+ fe 2)/4 a (4t) 



a(4i) a(4t)( fc i+ fc =)/ 4 9a; Va(4i)W4 
(fcx + fc 2 ) a (4t)/(^L_) J 

a (4i)(fc 1+ fc 2 )/4+l 

k 2 a{At)x /( 0(4*^2/4) 



r/(- 



a (4t)(fc 1+ fc 2 )/4 ' V a ( 4t )fc 2 /4 ' a ( 4t J 

■^(^(4*7^274 ) a(4t) 



fc 9 

) ; J ., , : — d(4t) -4 



a(4*) a(4i)( fe i+' £ 2)/4 a (4 f )fe2/4 



+ (fci + fc 2 ) 



a (4t)(fci+fc 2 )/4 a ( 4t ) 



0. 



(20) 
(21) 
(22) 

(23) 
(24) 



(25) 



The proof is completed. ■ 

On the other hand, in [5] and |23j , the qualitative properties of the Emden equation 

V ' a(s) K 

o(0) = a ^ 0, d(0) = ai 
where for k > 1 in [5] and [23] and k = | in [25] ; 

were studied. Therefore, the similar local existence of the Emden equations |T 

V ' a{s) K 

a(0) = a 7^ 0, d(0) = ai 
can be proved by the standard fixed point theorem [5] and [53]. To additionally show the blowup 

property of the time function a(s), the following lemmas are needed. 
Lemma 4 For the Emden equation Uty) . 

~ 4w = 
o(0) = a > 0, d(0) = ai, 



(26) 



(27) 
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the solution exists locally. 

After obtaining the nice structure of solutions (flT) , we just use the techniques of separation of 
variable ([13], [IS], [TB], [22], [23] and [25]), to prove the theorem: 

Proof of Theorem [T] From Lemma [3] it is clear to see our functions (|15[) fit well into the mass 
equation, (|3J)i, except for two boundary points. 

The second equation of 2-component Camassa-Holm equations ([3J2, becomes: 

u t - u xxt + Auu x - 3u x u xx - uu xxx + k 3 pp x (28) 
= u t + 4uu x + k 3 pp x (29) 



As the velocity u, in the solutions (|15|) is a linear flow: 



we have 



The equation ([29]) becomes 



aUt) , , 



u xx = 0. (31) 



u t + 4u x u + k 3 pp x (32) 

e (mn x + 4 m) ««, + fe ikfe* f *igS> ^ (33) 



dt\a(4t)J \a(4t)J a(4t) o(4t)( fc i+^)/ 4 U(4t)C=i+^)/ 4 

( MAt) q(4t) 2 \ /q(4*)\ <x(4*) /(^PSTt) gggM 1 

V o(4*) a(4<)V W 4 *)/ o(4*) a(4*)( fe i+ fe 2)/4 a (4i)(fei+fe2)/4 a ( 4i )fe 2 /4 V° J 

,a(4f) / ( a (4t)*2/4 )/( a Ut%'* ) ,„_,. 

a(4*) a(4t)^+^r 

^T-V + mHv)) (36) 



a(4t)^ + 3 -¥ \h 
for k 3 7^ 0, with the Emden equation: 

a(s) - t4t7t = 

V / 4a(s) K 



(37) 

a(0) = a / 0, a(0) = a x 
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by defining the variables s :— 4t, r\ := x/a{s) k2 / A and k = ^- + k-2 — 1. 
(l)for &3 = 0, we have £ = : 

/, s Pol 7 ?) /, \ a (4t) 

p(t,x) = — — 7-7-. — —j — T-j-7 , ufr.x) = . , ■ X (38) 

where po is an arbitrary C 1 function. 

Now, we can separate the partial differential equations into two ordinary differential equations. 
Then, we only need to solve for ^- < 0, 



/(0) = -a < 
or for 4- > 0, 

£»7 + /(»/)/(»») = o 

/(0)=a>0. 

The ordinary differential equations (|39|) or (|40| can be solved exactly as 



(39) 



(40) 



/w=fV^ 2+ U a) - (41) 



In fact, we have the self-similar solutions in details: 

(2) For fc 3 > and f > 0, or 

(3) For fc 3 < and £ < 0: 

=max( a(4f) f fc ( i ,? H ) fc2)/4 , 0), u(t,a;) = ggz 
«( s ) ~ W = °' a ^ = °° ^ °> a ^ = ai (42) 

With the assistance of Lemmas HI we may obtain the local existence of the above solutions. 
The proof is completed. ■ 

After we construct the solutions (fT5|) . the blowup or global behavior can be analyzed from the 
Emden equation: 

a(s) - -A* = 
^ ' ° (s) (43) 

a(0) = a ^0, a(0) = a x . 
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To additionally show the blowup or global property of the time function o(s), the following lemma 
is needed. For the particular case of ki = 1 and &2 = 1, we have 

«=y + *2-l = i<l. (44) 

For the case k = 1 or n = 1/3, the blowup or global results for the Emden equation are already 
shown in 2'A and [25] . For the case < k < 1, we can show by the energy method [15]. Due to 
the proof is very similar, we omit the details here to have the lemma: 



Lemma 5 For < K < 1, the Emden equation U2\) . 

«(*)-z(!f = 



(45) 

o(0) = ao > 0, d(0) = a%, 



(V £ < 0; fiere exists a finite time S , such that 



lim a(s) = 0. (46) 



f^J if £ > 0, the solution a(t) exists globally, such that 

lim a(s) = —oo. (47) 

s— )-+oo 

After obtaining the above lemma, it is clear to have the following result: 

Corollary 6 For < n < 1 arid ao > 0, we have: 

(1) &3 > and £ > 0, £/ie solution H5\) exists globally. 

(2) ^3 < and £ < 0, the solution M5\) blows up on a finite time T. 

For the other cases, the interested reader may determine easily by the classical energy method. 
In particular, for k > 1, the blowup or global behavior may be referred [5]. 

Remark 7 For < ana 1 £ < 0, the blowup solutions h!5\) collapse at the origin: 

lim p(t, 0) = +oo (48) 

t->T~ 
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with a finite time T; 

For ks > and £ > 0, the global behavior of the solution \15)) at the origin is: 

lim p(t,0) = 0. (49) 

t— >+oo 

Remark 8 The solutions U5\) are only C° functions, as the function f(rj) is discontinuous at the 
two boundary points, for a > 0: 

lim f( V ) ± 0. (50) 

Remark 9 Our analytical solutions could provide concrete examples for testing the validation and 
stabilities of numerical methods for the systems. Additionally, our special solutions can shed some 
light on the understanding of evolutionary pattern of the systems. 

Remark 10 We may calculate the mass of 

(1) for do > 0, the solutions il5\) with (k% > and £ > 0) or (k^ < and £ < 0): 



Mass = / pit, x)dx < +oo; (51) 

J — oo 

(2) for a < and (-1)" = — 1, the solutions fffljl with (k 3 > and £ < 0) or (fc 3 < and £ > 0) : 

Mass = +oo. (52) 



3 Blowup Phenomenon for > 

In this section, we show the blowup phenomenon by the standard method. The system can be 
converted to a quasi-linear evolution of hyperbolic type: 

Pt = -k 2 p x u - (h + k 2 )pu x 

(53) 

ut + u ■ u x + d x G * (|u 2 + !fp 2 ) = 

where the sign * denotes the spatial convolution, G(x) is the associated Green function of the 
operator (1 — d 2 )^ 1 . For the local well-possness and the blowup phenomenon with odd initial 
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values (po: u o, ), are discussed in [TU] recently. We may show the blowup by the particle trajectory 
method. Consider the following problem: 

q t = u(q,t), 0<t<T 

(54) 

q(x, 0) = x. 

we apply the characteristic curve method of hyperbolic partial differential equations, as the 
following theorem: 

Theorem 11 Suppose the velocity u is uniformly bounded at some point xq, such that u(t,xo) < 
M and u x (x Ol 0) < —y~^\M\and the regularity of the solution at xq is sufficient enough. The 
nontrivial solutions blow up before a finite time T. 

Proof. In general, we show that the p(t,x(t;x)) preserves its positive nature as the mass equation 
(|53]) i can be converted to be 

Pt + k 2 p x u = -(hi + k 2 )pu x (55) 

^7 + {kx + k 2 )p^- • u = (56) 
Dt ox 



D d t d . 
TT, = IT. + fc 2 [u— . (57) 



with the material derivative: 

/.)/ 01 "'V'o.r 

We integrate the equation ([56]): 

p(t, x) = po(x o (0, x Q )) exp ( -(fci + k 2 ) I V • u(t, x(t; 0, x Q ))dt ) > (58) 



for po(^o(0, a^o)) ^ 0> along the characteristic curve. 

To drive the argument for blowing up, we differentiate equation (|53[) 9 with respect to x: 

2 3 2 k 3 2 /3 2 k 3 2 \ 

u x t = -u x - uu xx + -u + —p -G*\~u + —p J. (59) 

Along the another characteristic curve with 

d d d 

dt := dt- U dx-> {60) 



we have for some point Xq: 



du x (x ,t) _ 2 3 2 . / 3 2 . h ., 



iU K(x ,t) + -u \x .t)-G* \^-u z + ^p 2 ) (x Q ,t) ii.h 

<-ul(x ,t) + hl 2 . (62) 
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The above Racci equation blows up before a finite time T, if the initial value 



u x (x o ,0) < 




(63) 



The proof is completed. ■ 
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